% clear all ; clc ; close all

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
a=1./[10 100 1000 10000];
n=1.5;
% m=1;
m=1-1./n;

psi=logspace(-2,8);

for i=1:length(a)
    S(i,:)=(1 ./ ( 1+(a(i)*psi).^n ) ).^m;
end

figure
semilogx(psi,S);
xlim([10^-2 10^8])
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
a=0.1;
% n=[0.5 1 2 4];
% m=1;

n=[1.25 1.5 1.75 2];
m=1-1./n;

psi=logspace(-2,8);

for i=1:length(n)
%     S(i,:)=(1 ./ ( 1+(a*psi).^n(i) ) ).^m;
    S(i,:)=(1 ./ ( 1+(a*psi).^n(i) ) ).^m(i);
end

figure
semilogx(psi,S);
xlim([10^-2 10^8])